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Abstract 

We study the effects of an information-theoretically motivated non- 
linear correction to the Wheeler-deWitt equation in the minisuper- 
space scheme for flat, k = 0, Friedmann-Robertson- Walker (FRW) 
universes. When the only matter is a cosmological constant, the non- 
linearity can provide a barrier that screens the original Big Bang, 
leading to the quantum creation of a universe through tunneling just 
as in the k = 1 case. When the matter is instead a free massless 
scalar field, the nonlinearity can again prevent a contracting classical 
universe from reaching zero size by creating a bounce. Our studies 
here are self-consistent to leading order in perturbation theory for the 
nonlinear effects. 



1 Introduction 

When extrapolated back to early times, classical cosmologies typically lead 
to singularities and a hence breakdown of predictability in our physical laws. 
Consider the Einstein-Hiibert action for a FRW universe, 



$ = I dtL= l -\ (!tXa :i 



-a 2 (f) 2 . s k 



(1) 
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where a(t) is the dimensionles^l scale factor, V{<j>) the potential energy of 
the scalar field <j)(t), N the lapse function, k = 0,±1, and we have taken 
h — c—1. For slowly varying scalar fields one may ignore the kinetic energy 
term and set 1/V = Oq, so 

S = JdtL = ^JdtN 

Varying with respect to N and choosing the N = 1 gauge gives us the 
Friedmann equation for an isotropic and homogeneous universe with the 
cosmological constant A = 3/a 2 , modeling inflationary sources in the early 
universe, 

d 2 + (k - a 2 /al) = . (3) 
The flat, k = 0, geometry has the expanding classical solution 

a = exp(t/ao) (4) 

which implies an arbitrarily small universe, a — > 0, at early times t — > — oo. 

Of course when the universe is small there is no reason to expect classical 
physics to be accurate. At such scales quantum effects should be important 
and one hopes that they cure classical singularities if and when they occur. 
In quantum cosmology studies [2], the Wheeler-deWitt (WDW) equation is 
a functional Schrodinger equation that provides one with a wavefunctional 
to compute quantum effects in the universe. In the minisuperspace scheme 
the quantisation is applied to restricted situations such as the FRW model 
leading to an ordinary Schrodinger equation: From ([2]) one gets the canonical 
momentum £> = §f = —da. Promoting this to an operator, p — > p = —i-§^ 
and applying the quantisation to the Friedmann equation ([3]) leads one to 
the WDW equation in minisuperspace [3], 

#0=0, (5) 

which is just a one-dimesional, < a < oo, time-independent Schrodinger 
equation for a particle of mass m = 1/2 moving in a potential V(a) = 
a?{k — a 2 /al), with units % = 1. 

2 The physical scale factor a p h ys = /3a, where [3 = *J -$y-, l p being the Planck length 
and V is the volume of the spatial hypersurface divided by a 3 pQ . 
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For a spherical geometry, k = 1, the potential barrier between a = and 
a = ao means that there is quantum tunneling [H [3], allowing a classical 
universe to emerge and exist at finite size (a > ao). However in this case 
the zero size, a = 0, universe is actually already avoided in the classical 
dynamics, Eq.(j3j, as the classical potential creates a bounce for a classical 
universe starting from a > and evolving backwards. 

The classical dynamics is different when k = 0, as seen in (J3J), because 
there is then no barrier to prevent a backward evolving universe from col- 
lapsing to zero size. Standard quantisation of the k = geometry does not 
help: When one chooses the outgoing wavefunction as a solution to (JSj), rep- 
resenting an expanding universe, one still encounters the a = possibility in 
the sense that the wavefunction allows a universe of arbitrarily small size to 
exist, see Section(2). 

It is generally believed that the functional WDW equation is only an 
approximate, perhaps semi-classical, description of quantum gravitational 
effects, partly because of ambiguities in the functional fomalism. For those 
and other reasons, various modified WDW equations have been investigated, 
obtained, for example, through a postulated non-commutativity [7] or by 
using ideas inspired by loop quantum gravity [5j [6] . In the latter, a discre- 
tised WDW equation emerges which is found to avoid classical singularities 
through a bounce. 

It is as yet unknown which, if any, of the suggested modifications to 
classical general relativity is an accurate description of likely new physics at 
the Planck scale. We therefore adopt in this paper the "maximum uncer- 
tainty (entropy) principle" [8], or information-theoretic perspective, to mod- 
ify the WDW equation and study the consequences. The philosophy of the 
information-theoretic approach is that one should minimise any bias when 
choosing probability distributions, while still satisfying relevant constraints 



The maximum uncertainty approach is well known in statistical me- 
chanics [8]. For example, if one is only provided with the mean energy 
E = J e(x) p(x) dx of a classical statistical system, but not the probability 
distribution p(x), then one may infer the probability distribution by max- 
imising the Gibbs-Shannon entropy Iqs = — I p{x) kip(x) dx under the con- 
traint of the given mean energy. That is, one varies Iqs — f3E with respect 
to p(x) with (3 the lagrange multiplier. This gives the canonical probability 
distribution p(x) oc exp(— /3e(x)). Of course one must motivate the entropy 
(uncertainty) functional that is adopted, for example through an axiomatic 
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construction as was done by Shannon [5]. 

The above approach is now used in many fields and it may also be used to 
motivate the usual Schrodinger equation 0, [10] as a one parameter extension 
of classical ensemble dynamics, arising through a simultaneous minimisation 
of the Fisher information measure. The Fisher measure, 



is the simplest information (or inverse uncertainty) measure which satisfies 
axioms suitable for the ensemble dynamics context [TUJ , just as the Gibbs- 
Shannon measure is the simplest measure suitable for classical statistical 
mechanics. Deformations of the Fisher measure will then give rise to gen- 
eralisations of the Schrodinger equation suitable for probing potential new 
physics at shorter distances. It turns out that such deformations lead natu- 
rally to nonlinear extensions of Schrodinger's equation [TT]. The information 
theory approach also allows for systematic improvement of the models as our 
knowledge progresses. 

As we saw above, in the minisuperspace scheme Eq.flS]) is just the time- 
independent Schrodinger equation defined on the half- line < a < oo. We 
postulate therefore that the unknown new physics at shorter distances may be 
modeled, within the information theory framework, by a nonlinear correction 
as found previously in Ref. [TT], so that the modified equation becomes 





da 2 



+ V(a) + F(p) V(a) = 



(7) 



where 



F(p) = Qinl-Q, 



(8) 



with 



QlNL 



1 
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+ 1 - 



(1 - rj)p 







(9) 



and 



Q 



1 d 2 ^/p 
yfp da 2 



(10) 
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Here p(a) = ip*(a)ijj(a) and p±{a) = p(a ± 77L). As we have made a di- 
mensionless, L > is the dimensionless (scaled by (3) nonlinearity scale and 
< r] < 1 is a parameter that labels a family of nonlinear isationqj . 

At the level of the action, F(p) is obtained by varying the Kullback- 
Liebler (KL) information measure Ikl which in the limit L — > reduces to 
the Fisher information measure responsible for the usual linear Schrodinger 
equation [IT] . Here is a brief review of the basic KL measure given by 

/p(x\ 
p(x)ln——dx. (11) 
r[x) 

It measures the relative uncertainty between two probability distributions, 
p(x) and a reference r(x). If we have no useful a priori information then r(x) 
can be taken to be a uniform distribution and the KL measure then reduces 
to the Gibbs-Shannon entropy. On the other hand if in Eq. ffTTl) one chooses 
the reference distribution r(x) to be the same as p(x) but with infinitesimally 
shifted arguments, that is r(x) = p(x + Ax), then to lowest order, 

-(Ax) 2 

Ikl(p(x), P (x + A(x)) = ^-^I F (p( x )) + 0{Axf. (12) 

So to lowest order, minimising the Fisher information, which leads to the 
linear Schrodinger equation, is the same as maximising the KL uncertainty 
measure for two probability distributions that are close to each other, r(x) = 
p(x + Ax); one might interpret Ax as the resolution at which the coordinates 
become distinguishable [11]. Thus the KL measure interpolates between the 
simplest inference measures used for quantum theory and statistical mechan- 
ics. 

Furthermore, the (regularised) Kullback-Liebler measure used to obtain 
([7j) is arguably the simplest generalised measure which has properties similar 
to those of the Fisher measure and yet leads to nonsingular equations of 
motion. In particular, the nonlinear correction F preserves the invariance 
of the Schrodinger equation to changes in the norm of the wavefunction, 
ip — > Xip. We refer the interested reader to [SHU] for further details. 

In summary, the motivation for choosing the regularised KL measure, 
and hence the family of nonlinear equations ([7]), is that it leads to a mini- 
mal deformation of the Schrodinger equation within the information theory 

3 As the equation ((7J is only defined on the half-line, we do not symmetrise the nonlin- 
earity with respect to L [llj . using the the simplest version here. 
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framework. As stated above, we intend to study the consequences of these 
nonlinear WDW equations, focusing in this paper on some simple situations. 
In the next section we will show that the nonlinearly corrected quantum 
dynamics can avoid the zero size, a = 0, possibility that is present in an 
evolving flat FRW universe with a cosmological constant, Eq.(j4]). Then in 
Section(3) we study a flat FRW universe where the matter is represented by a 
free massless scalar field. Here we show that the original classical singularity 
can be avoided by the modified classical equations induced by the nonlinear 
correction to the quantum dynamics. 

We remark that the issue of singularity avoidance in cosmology has been 
previously studied by several authors using different approaches, see for ex- 
ample the reviews [51 El EE2] arid references therein. What we are proposing 
here is a different framework, that of information theory, to motivate any 
corrections to the dynamics of the universe. As discussed in the concluding 
section, we hope to use this framework to also study other issues in cosmol- 
ogy and to see if this method of "minimum bias" can be related to other 
approaches in the literature. 



2 Cosmological Constant 



The flat, k — 0, universe is favoured by current data so it is of interest to see 
if the nonlinear correction helps screen the initial singularity typical for such 
FRW universes. Here we study the case when matter is only in the form of a 
cosmological constant. Setting k = and a = lb, with I = a^ 3 , in ([7]) we get 

b A + l 2 F(p(lb)) 



db 2 



0(6) = 0, 



(13) 



ijj(a) = 4>(b). We assume that the nonlinearity is small even in the early 
universe so that F may be expanded perturbatively in L to lowest order, 

77(3 -Arj)L 



F(b) 



where 



and q{b) 



f(b) 



I 3 
q' 



12q 3 

'(b)(j)(b). Thus we may write 
d 2 



f(b) + 0(L 2 ) 
(2q' 2 - 3q"q) 



db 2 



6 4 + 77(3-477)6/(6) 



0(6) = 



(14) 
(15) 

(16) 



6 



where e = L/l measures the strength of the nonlinearity. Assuming e < 1, 
we may solve the equation by iterating about the unperturbed solution for 
e = 0. The unperturbed solution which represents an expanding universe at 
large scale is given by a Hankel function 

Mb) oc VbH[%(b 3 /3) (17) 



nb 3 



exp 



(> i(6 3 - tt)/3 as 6^oo. (18) 



Recalling p = —ad, the negative momentum of this solution corresponds to 
a > as required. From this unpertubed solution one calculates go = 0o0 
and then fo(b) which is shown in Fig.(l). Then to leading nontrivial order in 
the iteration one has a linear Schrodinger equation with an effective potential 

V eff = -6 4 + 77(3 - 4*7)6/0(6) • (19) 

For small b, 

fo « 0.16 (20) 

and so for rj < 3/4 there is an effective potential barrier, a finite size universe 
coming into being through quantum tunneling. 

It is interesting to calculate the probability for tunneling through the 
barrier. For the k = 1 case the result in the WKB approximation is known 
in the linear theory [3J, 

P k=1 = exp [-2a^/3] , (21) 

with small nonlinear corrections not affecting the result significantly, so the 
universe is likely to be born with a small size [3J . For the present k = case 
the WKB formula applied to the effective Schrodinger equation gives 



J fe=0 w exp 



dby/V eff (b) 



(22) 



with bi the point where V e ff(bi) = 0. Using the expression (120!) we may 
estimate 

61 = 0.46(r7(3-47/)e) 1/3 (23) 

and hence 

P k=0 ^ exp [-0.177(3 - 477)e)] , (24) 
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so that for fixed rj < 3/4 one may say that small values of e are "preferred" , 
a conclusion which is self-consistent with our approximation e 10. 

Having obtained the effective potential above, one may alternatively dis- 
cuss the modified classical dynamics of the universe. It is clear now that 
for parameter values rj < 3/4, a backward evolving classical universe will 
experience a bounce instead of shrinking to zero size. 



3 Free Massless Scalar Field 

As the WDW equation is independent of time, so in quantum cosmology one 
often looks at correlations between variables to describe the evolution of the 
universe [2] . One possibility is to use a free massless scalar field as an internal 
clock [6J, which we also do below. The coupling of a free massless scalar field 
to gravity may also be viewed as the other extreme to the model we studied 
in the last section; here we have kinetic energy but no potential energy for 
the matter. 

Following convention^, we define a = In a so that the k = classical FRW 
action becomes 



S i dtL= X -\ ,11 X< ;: " 



-a 2 



N 2 + N 2 



(25) 



giving the following classical equations on varying with respect to a,<f>,N 
and then setting N = 1, 

+ 30d = 0, (26) 

2a + 3d 2 + 30 2 = 0, (27) 

-a 2 + <p 2 = 0, (28) 

the last equation being a constraint. After a scaling of the time variable, we 
have a solution 

a = ^hgt + C, (29) 
= ^logt (30) 



4 The k = — 1 results will be discussed elsewhere. 

5 Actually we need this transformation to get the Klien-Gordon equation ([32]) with con- 
ventional kinetic terms for both variables so that we may implement the nonlinearisation 
as in Ref . 1111 . 



8 



where C is an arbitrary constant (we have only indicated the classical solu- 
tion which interests us, for which a increases as the internal time <fi increases). 
As t — > 0, a — > — oo and so a — > 0, the classical universe arises from an initial 
singularity. 

To quantise the classical theory, we start as usual from the Hamiltonian 
corresponding to 



H 



Ne 



-3a 



2 , 2 
-Pa + P<p 



(31) 



where the canoical momenta are p a 



'* and 



x r<p = ^-jf-- The classical 

constraint (|28j) becomes the Hamiltonian constraint Hip = giving us the 
WDW equation 

d 2 d 2 ' 



+ 



da 2 



ip(4>, a) — . 



(32) 



The general solution of (132]) can easily be obtained by separation of variables 
and then a wavepacket constructed [T3], [I], 



w{k) A k {a)B k (4>) 



where w(k) is an arbitrary function of k and 

Ak{a) = a\ exp(ika) + a 2 exp(— ika) , 
Bk{<P) = bi exp(ik(p) + 6 2 exp(—ik(p) , 



(33) 



(34) 
(35) 



for some constants ai,a 2 ,bi,b 2 - We would like the state to represent a 
large universe when the intrinsic time <ft is large, and so we construct such 
a localised wavepacket by taking a 2 = b\ = and the Gaussian weight 
w{k) = exp — ( k ~!p to get 



ip = aib 2 <j\/7r exp 
a and d being constants. The density 

p((j),a) = ip*4> = (ai620") 2 7rexp 



exj>(id(a — 4>)) 



(36) 



a — 



(37) 
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is clearly localised near a ~ for large 0; we take this to represent a classical 
universe at large time. But we also see that the universe would have been 
arbitrarily small, a — > — oo, in the distant past <ft — > — oo. 

We now proceed to nonlinearise the Klien-Gordon equation (132]) following 
the information-theoretic approach discussed in Ref.JTT] to get 



d 2 



+ 



d 2 



d(f) 2 da 



-F a (p) + F^(p) 



ip(<p, a) = 



(38) 



where F a and have the same form as the F in dSJ) but with generally 
distinct nonlinear parameters L Q > and > corresponding to the 
gravitational and matter degrees of freedom. 

As in the previous section, we solve (1581 by perturbation and iteration to 
lowest nontrivial order using (j3~T|) to get an effective linear equation 



with 



d 2 



+ 



d 2 



d<p 2 da 



eff 



a) = 



(39) 



V. 



eff 



U 



U 



a' 



a) 2 -3 



— a) 



7/(3- 477) (L a + L^)(T 4 
12 



(40) 
(41) 



This linear equation describes the approximate, to leading nontrivial order 
in perturbation theory, quantum dynamics of the states that are still highly 
localised. 

As the hyperbolic nature of the Klien-Gordon equation makes the tun- 
neling analogy problematic, see for example Ref . [14] , we proceed differently. 
We obtain the effective classical equations that imply, through the correspon- 
dence principle, the modified quantum equation (|39|) . First we note that the 
modified equation (1591 corresponds to the Hamiltonian 



H 



Ne 



-3a 



-fa+Pl + V eff(<P,a) 



(42) 



and this arises from canonically quantising the classical action 



,s : ■ / dtL= l -\ <ltXr'"> 



—a 



iv 2 + h~ e ~ 6av ^ 



(43) 
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This effective classical action gives the following modified evolution equations 
(in the N = 1 gauge) 



Note that 





6 + 36a + -e ^ = 


o, 


(44) 


2a + 3d 2 + 3 


2 + e~ 6a 


Uv dVeff ~ 

Weff da 






(45) 




-a 2 


+ <j) 2 + e~ 6a V eff = 


0. 


(46) 


dVeff _ 

86 


dVeff _ 

da 


3u a 2 {6 — a) 2 


- 1 




(47) 



It is easy to verify that the modified constraint equation (|46p combined 
with any of the other two evolution equations (I44II45P implies the third. With- 
out loss of generality in the subsequent analysis we set a = 1. 

We remark that since the quantum states used to obtain the effective 
potential fl4ll were localised, the effective classical equations (I44H46I) are a 
self-consistent description of the mean dynamics of such states. 



3.1 Analytical Results 

A subset of solutions to the new coupled equations (1441146 p can be studied 
analytically by assuming the correlation a = 6 at all times. The constraint 
equation is then automatically satisfied and the other equations reduce to 

3a 2 + a- —e- 6a = (48) 

which is easily integrated to give 

e 6a a 2 = 3ua + C (49) 

where C is the integration constant. Since the left-hand-side is positive, 
therefore taking rj < 3/4 as before (that is, u > 0) gives 

—C 

a> (50) 

~ 3u K ' 

and so 

— C 

a > a min = exp(— ) (51) 
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That is, there exists a class of solutions for which a(t) = <f>(t) for all t: 
Such universes have a minimum nonzero size, which depends on the initial 
conditions that fix C, and hence the Big Bang singularity is avoided, being 
replaced by a bounce. Note that while a a = <fi correlation is also possible 
in the original classical dynamics, (I29|30p . only in the modified dynamics is 
the singularity avoided. 

The a = <ft solutions of the modified dynamics are stable to small pertur- 
bations as one can see by setting a(t) = <p+5(t) and linearising the dynamical 
equations for 5 1 to get 

* + 2sSr° (52) 

which is the equation for an oscillator with time-dependent frequency; 5 is 
thus bounded. 

It should be emphasised that though the a ~ <fi solutions to (14411461) form 
a subset of solutions to those equations, they are a very important class 
because the effective potential (pfOj) was constructed precisely from quantum 
states (1371) which were localised around a ~ 0. Furthermore, to reiterate, 
while in the original dynamics such states encounter a singularity, in the 
modified dynamics one sees instead a bounce. 



3.2 Numerical Solution I 

We now investigate the coupled effective classical equations ( |44T|4T)|) numer- 
ically. We show in the Appendix that if one solves the two second-order 
equations (144ti45j) with initial conditions at t — to that satisfy the constraint 
equation (146!) . then the constraint equation is satisfied for all times. The 
initial configuration is chosen to be a large universe (a > 1) at large time 
t. Notice that the nonlinear corrections come with the pre-factor exp(— 6a) 
which causes the nonlinear terms to be negligible when the universe is large. 
So our initial configuration satisfies very closely the classical, unmodified, 
FRW equations which we had considered earlier fl29|30l) . 

Consider first the choices, of initial time to — 27 and C — 0. This gives 
initial conditions «o = 0o = log 3 and do = 1/81, with 0(to) determined by 
the constraint equation. We set u = 0.01. The numerically solved equations 
give the plot shown in Fig. (2). The solution is found to satisfy a = <fi for all 
time with a min ~ —2.5. In the figure we also compare the modified trajectory 
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with the classical result, showing a bounce in the first case and a singularity 
in the latter. 

In the Appendix we compare this numerical solution with the analytical 
results from the previous subsection. 

3.3 Numerical Study II 

By choosing nonzero values of C in (|29|) we can study more general initial 
conditions for which a ^ f We keep to = 27 and u = 0.01. Figs. (3-6) 
show some examples, in all cases a bounce occurs at some nonzero size that 
depends on the initial conditions. 

As we keep increasing a beyond log(3.5), while keeping t = 27 fixed, the 
curve for a(t) shows an increasingly sharper turning point and increasingly 
lower minimum. This might seem to indicate that a singularity, a — > — oo, 
will eventually appear when starting from some initial conditions correspond- 
ing to large for a fixed t$. However note that in those cases where a be- 
comes increasingly negative, the corrections from V e ff to the classical equa- 
tions have become large, recall the prefactor exp(— 6a), and we are thus 
beyond the validity of our perturbative approximation. Therefore all that 
one can conclude at this point is that within the regime where the perturba- 
tive approximation is valid, the singularities are resolved, being replaced by 
a bounce. 

4 Conclusion 

We have used an information-theoretic approach to model the unknown new 
physics that is generally believed to exist near the Planck scale. For simplic- 
ity, the resulting nonlinear WDW equation in minisuperspace was analysed 
assuming the nonlinearity to be weak. The nonlinearity in the WDW equa- 
tion induces an effective potential in the modified classical dynamics: In the 
modified k = FRW models with a cosmological constant or a free mass- 
less scalar field, a classically contracting universe may experience a bounce 
instead of reaching zero size. 

The nonlinearisation of the quantum equation makes the Hamiltonian, 
and hence the effective classical potential, state dependent. The state- 
dependence is not necessarily undesirable: It could perhaps be used to se- 
lect states that lead to mathematically well-defined evolutions and hence 
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hopefully guide us to the effective (emergent) laws that govern reasonable 
universes like the one we inhabit. 

The state-dependence of our bounce is in contrast to some other ap- 
proaches such as that of loop quantum cosmology [6|, [5] . In passing we note 
that in loop quantum cosmology the modified WDW equation is a linear dif- 
ference equation while in our approach it is a nonlinear difference-differential 
equation; this seems intriguing. 

In the free massless scalar field case, the effective potential that was gen- 
erated through the nonlinearity was responsible for the bounce and hence 
also a momentary acceleration of the universe. However as we used a pertur- 
bative approach, the inflation is mild. It would be of interest to study the 
nonlinearity non-perturbatively to see if realistic inflation can be generated 
by the nonlinearity starting from free matter fields, as there are indications 
that the nonlinearity has novel non-perturbative properties |15j . Indeed, one 
might expect the nonlinearity to be large near the Planck scale and perhaps 
weaken as the universe expands, so it might be useful to make the nonlinear- 
ity parameter L dependent on the scale factor a(t). 

As we studied the realisation of the bounce in two extreme cases of matter, 
pure potential energy and pure kinetic energy, we think that similar conclu- 
sions will hold for more general matter though this would require checking. 
It would be useful to study in a realistic cosmological model not only the 
singularity avoidance but also the residual long time consequences of the 
nonlinearity, for example its effect on the cosmic microwave background and 
its relevance to the present acceleration of the universe. Such phenomenol- 
ogy might help fix, and perhaps also provide physical meaning to, the free 
parameter < i] < 1 in our model: Recall that only for rj < 3/4 did we 
find a bounce in the perturbative regime. Phenomenology would eventually 
also be important for discriminating between the predictions of the various 
different approaches to quantum cosmology in the literature. 

Of course one should go beyond the FRW minisuperspace and study less 
symmetric geometries that are better models for a beginning universe. Other 
gravitational systems for which an information-theoretically modified dynam- 
ics might yield interesting results are quantum black/worm holes and the 
current classical accelerating universe. 
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Appendix: Numerics 
A Lemma 

Denote the left-hand-side of the constraint (j46l) by 

A(t) = -a 2 + <p 2 + V(a, (f)) (53) 

where we have also defined V(a, <fi) = e~ 6a V e ff. Multiplying fj44"l) by <f> and 
f )45|) by a and using those results in A, one obtains 

A = -3aA (54) 

or 

e 3a A = constant. (55) 

Thus, when solving ( 1441451) simultaneously, if the initial conditions are chosen 
such that A = at some time to, then the constraint equation ( H6l) is satisfied 
for all time. 

The consistency of the constraint equation may be discussed more gener- 
ally using constraint algebra arguments, as for example in Ref.|16j. 

First Order Constraints 

The numerical solution of (|44"l 145]) was obtained using the Maple 11 software. 
Initially we used the Fehlberg fourth-fifth order Runge-Kutta method and 
rechecked the results using the forward Euler method. 

In addition we performed the following cross-check to detect numerical 
errors. The two equations (I44|45p are two second-order equations. By com- 
bining them in certain ways we obtain a new first-order constraint as follows: 
Multiply Eq.flU by 3/2, add Eq. (j24l . subtract Eq. fpl) and use Eq.(|4T]) to 
get 

x + 3xa = (56) 



15 



where x = (f) — a. Integrating the last equation gives 



e 3a (0 — a) — constant. 



(57) 



Together with the original constraint equation (1461) . this new first-order con- 
straint can be used for checking the numerical results: The numerical solu- 
tions of f)44|45l) give us the values of a, <fr and their first derivatives at any 
point. Such data from several points in our plots were inserted into the two 
first-order constraint equations to check for consistency. 

Checking a a = (j) solution 

The solution discussed in Section(3.2) was found to have a = <p for all time. 
It had a m i n ~ —2.5 and t min fa —1.4. These values can be compared with the 
analytical results of Section(3.1). The constant C in (1491) can be evaluated 
using the data from Section(3.2) to give C ~ 0.078 and hence (15D|) gives 
otmin ~ —2.6 in good agreement with the results from the plot Fig. (2). Next, 
one may integrate fj49|) . 



Using the known parameters, this yields t min fa — 1.4 again in good agreement 
with the plot Fig. (2). 
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Figure Captions 

• Figure 1 : The induced potential fo(b) in Eq.(19) is always positive; it 
vanishes at the origin, peaks around b = 1.2 and approaches zero as 
the dimensionless scale factor b approaches infinity. 

• Figure 2 : For the initial conditions described in Section(3.2), a = 
loga = log(3): The dashed line represents the trajectory of the classi- 
cal k = FRW universe with a massless scalar field; it hits the singular- 
ity at t — 0. The solid line is the trajectory in the modified dynamics, 
showing a bounce. The curve for coincides that for a. (Natural units 
are used as mentioned in Section 1). 

• Figure 3 : Results for «o = log(2.5): The dashed line represents the tra- 
jectory of the classical while the solid line corresponds to the modified 
dynamics. 

• Figure 4 : Results for ao = log(2.5): Trajectories for a and in the 
modified dynamics. The curve is above the a curve at large positive 
time but below it for large negative time. 

• Figure 5 : Results for a = log(1.5): The dashed line represents the tra- 
jectory of the classical while the solid line corresponds to the modified 
dynamics. 

• Figure 6 : Results for ao = log(1.5): Trajectories for a and in the 
modified dynamics. The curve is above the a curve at large positive 
time but below it for large negative time. 
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